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1. Introduction 



In this paper, we present results of studies of efficiency of tfie critical- component 
direct method proposed in [1 3] for solving degenerate and ill-posed systems of linear 
algebraic equations 

(1.1) AZ^F, 

where A is a square matrix of the general form with real elements a^j,A = {a^j},Z 
is an unknown vector with coordinates Zp Z = {zj}, and F is a known vector with 
coordinates f^,F = {fi},i,j = l,2,...,m. It is shown that for systems like (1.1) the 
critical-component method makes it possible to numerically determine the only normal 
pseudosolution (Z+= A+F) : \\AZ+- F\\ = mf\\AZ - F\\,\\Z^ \ = inf ||Z||, where Z^ 

is a set of all pseudosolutions to system (1.1), and to obtain the unique matrix A'^, 
pseudoinverse of A: \\A+A-E\\^ inf || - F||, ||/l+|| = inf \\ A~^\\,A+A^ AA+, 

o 

where is a unit matrix and fl^ is a set of all A , pseudoinverse of A. In this case, 
even if the problem (1.1) is substantially ill- posed, the quantities Z~^ and A'^ are stable 
to small changes of input data {A,F). Comparative analysis of results of the numerical 
solution performed for a large number of problems like (1.1) both by the new method 
and by those known earlier shows that the critical-component method is on the average 
more effective than any method compared to it. When detA ^ and a system is well- 
posed, the normal pseudosolution Z"*" of system (1.1) derived by the critical-component 
method coincides with its usual solution Z, and = A^^ is a matrix inverse of A. One 
of the main problems in numerical solution of ill-posed systems of algebraic equations 
is well-known [4,5,6]: there can be large changes in the solution, beyond the scope of 
admissible values, corresponding to small changes in the matrix of a system or/and its 
right-hand side. The above breakdown of continuity of the inverse mapping Z = A~^F, 
if A^^ exists, is caused by a great norm and, as a result, by large ^ = cond A, 

the condition number of the system matrix (/x = ||^|| • ||/1~^||, if detA ^ and fi = oo, 
if detA = 0, where || ■ || are the corresponding norms), i.e. even for an exactly given 
vector F a negligible relative error in calculating A~^ can produce a large distortion of 
the searched vector Z. This effect is to be taken into account since realistic calculations 
are carried out with a certain finite accuracy and , besides, sometimes one knows not 
the exact system AZ = F, but only a system AZ = F, approximate of it, which obeys 
the inequalities ||Al — A|| < h* and ||F — F|| < S* (the meaning of norms is defined by 
the character of a problem). The numbers h* > and S* > 0, specifying the norms 
of deviations of approximate data {A. F) from the exact ones {A,F) of problem (1.1) 
{h* < ho + hi, 6* < So + Si, ho > 0,hi > 0,So > 0,Si > 0), are sums of (/iq, ^o), proper 
model (complete) errors of problem (1.1) and of {hi, Si), round-off errors [7,8] when writing 
the data into the computer memory. Since there are, thus, infinitely many systems (1.1) 
with the input data {A,F), indistinguishable within the accuracy {h*,6*), we can speak 
only about deriving an approximate solution to system (1.1). As a result, difficulties may 
arise in numerical computations for some systems of equations (1.1) with square matrices 
when answering the following questions: 
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- is the system degenerate "within accuracy {h*,6*)" ill-posed?*^ and 

- is a given system ill-posed by virtue of its being degenerate or is it nondegen- 
erate but ill-posed? 

Indeed, if the system AZ = F with a square matrix is degenerate, then dct A = 0, 
i.e., the matrix A has some of its eigenvalues equal to zero. But if detA ^ 0, and the 
system is ill-posed, then the normal matrix A has some eigenvalues only close to zero 
fif, fifn are singular values of the matrix A). Consequently, systems of linear alge- 
braic equations with square matrices, which are ill-posed and degenerate "within a given 
accuracy {h*,6*)" may turn out to be indistinguishable in the process of computations. 
Besides, the problems (1.1) and AZ = F can be inconsistent if one defines the criterion 
of consistency [9] determined by accuracies {h*,5*). It may also happen that detA = 
(or detA = 0), i.e. system (1.1) (or AZ — F) has an infinite number of solutions. Then, 
there arises the question: what is to be understood by the numerical solution to the initial 
system AZ = F. There are various conceptual approaches to solve this problem (see, for 
instance, reviews given in [4,6,10], etc.). 

If one takes advantage of the regularization [4] , the solution Z'^ to the system AZ — F 

(1.1) will be the regularized normal pseudosolution that minimizes the discrepancy 

I \AZ — F\ \ on the set of all its pseudosolutions if I I = inf \\Z\ \ and is stable to 

zez/ 

small variations in [h*, 5*) of input data {A, F). The parametric vector is directly com- 
puted by solving the sequence of normal systems of equations {A^A + aE)Z°' = A^F with 
the aim of a more accurate iterative determination of the minimum of quadratic functional 
M°'[Z,F,A] = \ \AZ — F|p + ci;||Z|p with the regularization parameter a{a > 0), deter- 
mined from the discrepancy, i.e., from the condition \\AZ°' — F\\ — 5^, where 5* (5* > 0) 
is a numerical function of (/i*, 6*) and Z" [4,5,6]. 

The other group of numerical methods of solving the problem (1.1) rely on searching 
for the generalized matrix A"*", which is (pseudo)inverse of A, either by the method of 
singular decomposition (A — C/EV, where U and V are orthogonal matrices, E is a 
diagonal matrix, whose elements are singular numbers > I/U2I > ■ ■ ■ > > of the 
matrix A, and A+ = V'^H'^U'^), or by some other method [7,9,10,11]. Common to both 
of the approaches is that in their program realization they solve (each by its own means 
and with its own efficiency) the problems of minimization of norms \ \AZ — F\\ and ||Z|| 
and of the continuous dependence of the solution Z^ on small changes in {h*,S*) of input 
data {A, F). Here it is set that fi = cond A = jjAj] • ||A"^|], and the main problem now is 
a stable calculation of the rank of A [7,9]. 



*^Systems degenerate "within accuracy {h*,6*)" are not always ill-posed (as example is system (1.1) 
with A = A'^, singular (eigen)values /Ui = /U2 = ••• = Mm = 10~^, determinant detA = 10~^™ and the 
condition number cond A = = • = 1). 
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Conceptually, the critical-component method can be attributed to the second of in- 
dicated groups of methods. It is based on the idea of constructive search (under the 
condition that matrix and vector norms are consistent: ||2'+|| < \\A^\ \ ■ \ \F\\; and the 
matrix norm is induced by the vector norm: ||A+|| = sup (| jA+Fl |/| |F| |) [5,9]) of an op- 

timal representation for the matrix pseudoinverse of the matrix A, in the process of 
decomposition of system (1.1) into subsystems, whose solution is stable to errors (^i, eo) 
and small {h*, S*) changes of input data (^4, F). High efficiency of the critical-component 
method is provided by its basic constituents: 

- the reduction, stable to errors {h*,S*;ei,eo), of system (1.1) to two- (tri) diagonal 

systems; 

- generalized processes {A, G}, stable to errors (^i, Eq) [14], for calculating ratios of 
upper (lower) corner minors of triangular matrices which allow one, accurate within 
constants {si and Eq) of the computer arithmetic, to determine the structure and 
diagonal elements of matrices that are inverse of them (introduced in [12,13]); 

- the algorithm of optimal (with (si, Eq)) decomposition of the system AZ = F into 
well-posed subsystems; 

- the algorithm of optimal sewing of the solution Z+ to the system AZ — F from 
well-posed subspace solutions. 

In what follows, along with problem (1.1) of the general form, we will consider the 
problems of numerical solution of degenerate and ill-posed systems of linear algebraic 
equations 



(1.2) 
(1.3) 



with square real matrices C3 and C2 of order m, of the tridiagonal and two-diagonal form 
respectively: 



(1.4) 



C, 



P2 Q2 



Pm—l^m—l^m 
Pm Qm 



<l2 



where X = {xi,X2, ...,Xm)'^ and X = 
(1/1,2/2, •••,2/m)^ and Y = (yi, ^2, ^r, 



[Xi,X2, 



are unknown vectors, and Y 



are given m-dimensional vectors, {gj}™^ are 



diagonal elements and {pi,ri}^2 are sub (off) diagonal elements of matrices C3 and C2, 

*•* Throughout we use the notation: ei{ei > 0) is the modulus of relative error of the arithmetic of 
computer operations with real numbers with a floating point; eo(£o > 0) is the modulus of absolute error 
of the computer zero 9, i.e. of any small real number (except for 0) from the interval G (0 — £0, £0 + 0), 
where is the usual zcroth element of the real axis, li 9 E (0 — £o,£o + 0) and ^? 7^ 0, it is accepted 
that 6 = [8,9]. Using constants £1 and £0, one can estimate [7] errors of arrangement (writing) of the 
real \m, m] matrix A and m-dimensional vector F in the computer memory in the form | \Acomp — A\\e< 



(ei\\A\\E+som = hi),\\F, 



comp 



■ F\\e< (£i||F||£; + soy/m = 5i), where || • ||b are the Euclidean norms of 



matrices and vectors, and hi > 0,Si> 0. 
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respectively. Without loss of generality, we consider systems (1.3) only with the right 

two-diagonal matrix. 

Since these problems are a particular case of problem (1.1), all said above applies also 
to problems (1.2) and (1.3), whose solutions and X+ are constructed more easily than 
Therefore, in the course of program realization of the above conceptions of solution 
of problem (1.1), the initial stage [4,5,7,9,13,15] consists in its reduction to problems (1.2) 
and (1.3), i.e. 




C2{Q^Z) = PF, where PAQ is a two-diagonal matrix, HA^A^, 
C3{Q'^Z)^Q'^F, where Cs^ Q'^AQ is a tridiagonal matrix, if ^ = A^. 



Here U'^U = E = UU^ , U: Q,P are matrices of reflections or rotations. The orthogonal 
transformations (1.5) stable*^ to errors {h*,d*;si,So), do not often improve the nature 
of the problem being ill-/well- posed. Ill-posed systems of type (1.1) sometimes can 
numerically be reduced to ill-posed systems of type (1.2) and (1.3), with the notation in 

(1.5) : X = Q^Z,Y = Q^F and X = Q^Z,Y = PF. Therefore the basic problem is 
numerical solution of such degenerate and ill-posed systems. Once the vectors X and X 
are obtained, we determine the solution to system (1.1), vector Z, in the form 

(1.6) Z = and Z = QX. 

Numerical solution of ill-posed systems (1.2), (1.3) with tridiagonal and upper two-diagonal 
matrices can be best reahzed by the following methods [4,7,9] : the inverse substitution 
with normalization, regularization, a singular decomposition with exhaustion. In sect. 3, 
we present (in particular) the results of comparison between computations performed by 
these methodes and by the new one. 



*)lt is known [5,7,9] that the Euchdean and spectral norms of matrices are invariant (theoretically) 
under the orthogonal transformations (1.5), i.e. there hold the equalities: ||C2||s = ||P^Q||_e = 
\\A\\e,\\C3\\e = WQ^AQWe = \\A\\e; IIC2II2 = \\PAQ\\2 = H^lh, ||C^|2 = WQ^AQW^ = \\A\\2 and 
II-PII2 = IIQII2 = 1, \ \P\\e = \\Q\\e = Vrn- As a result, n = cond A = cond C3 (or cond C2). Here 
||-P-F||e= ||y||_E, HQ'^^IIe = II^IIb; ll^lh is a norm induced by the Euclidean vector norms \\Z\\e and 
M{A) = m- max \aij\ and \\A\\e are norms consistent with norms \\Z\\e and However, 

l<2,j<m 

in real computations in process (1.5) of the reduction of system (1.1) to form (1.2) or (1.3), using the 
Householder U transformations (reflections), we obtain the estimates [7] 

\\{C2)com,- C^\e = WiPAQUmp- PAQ\\e < {[ (/_Tj-tl PI |i. + ^T^^^^^t ) ^ M^^W^ ^^4' 

WiPFUmp- PF\\e< {er\\F\\E + 0r)=52 and 

\\iC3Un^p-C^\E^\\iQ^AQU^,-Q^AQ\\<{[[J^^^ 

\\{Q^F)eomp-Q^F\\E < (£r||J^||B + 0^) = 62, whcrc Er ~ 29£i and (2m + 2V^)£o,/i2 >0,S2> 0. 

Similar inequalities could also be written for A, F, where matrix A and vector F differ from A and F 
by simultaneous inclusion of inherited errors and errors of writing into the computer memory. From the 
above inequalities it follows that problems AZ = F and AZ = F are continuous with respect to the 
orthogonal transformations (1.5). Though the inherited errors (ho,So), if known, are, as a rule, much 
larger than the total {hi + /12, Si + 62) effect of the errors of writing and transformations (1.5), the latter 
can influence the character (degree) of problem (1.2) or (1.3) being well-/ill-posed. The cited monographs 
contain also simplifled estimates for errors /12 and 82- 
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2. Critical-component method for numerical solution of degenerate and ill- 
posed systems of linecir algebraic equations with tri- and two-diagonal 
matrices 

Below we formulate the theorem according to which one can numerically obtain the 
only stable non-iterated normal pseudosolution X+ of the system of linear algebraic equa- 
tions of the general form (1.2), stable to errors (^i, Eq) and {h, 5), by the critical-component 
method. 

The vector X+ and the representation for the matrix consistent with it (Cg" = B), 

o 

pseudoinverse to C3, are determined as functions of stably computed vector X {& regular 

o 

component of X'^) and matrix B (a regular component of C3 ). In contradistinction to 
the problem of computation of singular numbers of matrices C3 being unstable in nature, 
the critical-component method is stable owing to the stable processes of computation of 
the ratios of upper (lower) corner minors {A, G} of this matrix. Thus, the method of 
solution based on the search for a non-parametric stable component of the pseudoinverse 
matrix [7,9] found one more argument for its being efficient (contrary to conclusions of 
perturbation theory according to which X+ and Cg" are not valid for computer calcula- 
tions) . 

Theorem. Let C^X = Y he either a degenerate or an ill-posed system of linear 
algebraic equations with a square, of order m, real tridiagonal matrix of the general form 
C3 (1.4). Also, let the system C3X — Y, where 11(73 — C3II < h and — y|| < 5, being an 
image of the system C3X = Y in the computer memory, be ill-posed but nondegenerate. 
Then the only pseudosolution X of the system C^X = Y that is minimal in norm (| \X~^\ \ = 
min), obeys the condition of the norm of discrepancy being minimal (| ICgX^— F 1 1 = min), 
and is stable to computation errors {ei, 60) and to small changes {h, 6) of the input data 
(C3, Y), can numerically be obtained by the following direct critical-component method*^: 

Start of computations: 
k — l,i — m; 

(2.1) k = i; 

[k] h [k] [k] i 0, if /c = 1, 

(2.2) x,= E B,^ y^, 0,= I [k] + , 

if i = /fc, then (2.5), otherwise (2.3); 

[k] 

(2.3) if I J < l/£i, then (2.4), otherwise A; = A; + 1 and (2.1); 

[fel 

(2.4) j = i + l,x,^+i=0; 

$ _ J \yj\- \Pj ^j-i '^j'+i I' at < 1, 

^' " I 1 I w , [^1 , [^1 I /I I . I I 1 



*^Here h< ho + hi + h2 and 6 < 5o + 61 + 62 if the system C3X = y is a reduced image of the system 
AZ = F; and h < Jiq + hi,6 < 60 + 61, where (ho > 0,6o'> 0) are hereditary errors and (hi > 0,6i > 0) 
are errors of writing the system C3X = Y into the computer memory if system (1.1) is initially of form 

(1-2)- . _ . 

Since numerical solution is derived for the system C3X = Y that is, within accuracy {h,6), indistin- 
guishable from the system C3X = Y, for simplicity of the notation, the very algorithm of numerical 
method and its proof are given in the notation of the system C3X = Y, i.e., without if this 

does not cause misunderstanding. The requirement det 6*3 ^ of the theorem will be removed later. 
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if < 2£i, then (2.5), otherwise A; = A; + 1 and (2.1); 

, [k] [1^] 
(2.5) xt=x,+ (j),; 

ii i — 1, computations are over, otherwise i — i — 1 and (2.2); 

End of computations. 

Here: 

[k] [k] 

Bij (4+1 <: i < h, 1 ^ J ^ 4 and k = 1, 2, n) are elements of submatrices B of the 

O O o 

matrix B=C^ that is inverse of a well-posed matrix C3 of the form: 



(2.6) 

o 

C3= 



?1 ^^2 
P2 Q2 ^3 



_pi(Wi+i=i) 





%+l^/3+2 
i^Z3+29Z3+2^i3+3 



^^2+1 



[2] 

--C 



13+1 



ih—rn) 





Pi2+2^'2+2'"'2+3 
Pm 



'n+l+l=l 7 




«3 + l 



«2 + l 



Zl=m 



[2] 
5 



^2 f h=Tn^ 



[1] 



where 
(2.7) 



[fc+i] 



— Pik+i+i B ik+iik+fi^:+^+i, k -1,2, ..,n 1 



[fc+i] 



[k+i] 



and B i^+iik+i diagonal elements of submatrices B which coincide with 

the last diagonal elements of submatrices, inverse of wcU-posed submatrices C k+i 

separated by the method; and n is the number of separated subspaces. 
[k] 

Elements B^j are calculated [1] by the formulae: 
(2.8) 

n P^: if 1 ^ j < h lk+1 + l<i<lk, 

for all i from j<i<lk, if Aj = 0, for any j from Ik+i + 2<j <lk, 

[k] 

Oforall j from l<j<i, if Gj=0, for any i from Zfc+i-|-l<i<Zfe — 1, 



[k] 



n /^f, if /fe+i + 1 <i < j <ik, 



[k] 



for all i from Zfc+i + 1 < i < j, if Gj— 0, 
for all j from i < j < Ik, if Aj = 0. 

[fe] [k] 

Diagonal elements Bu of submatrices B and quantities a;^ in (2.8) are calculated [1] by 
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the formulae: 



(2.9) 





(Ai+1 


Bu- 


lk] 


[k] 


Bu- 


0, Bi- 


lk] 


lk] 


k Bii — 


0:Bi- 



Ik] 



Ik] Ik] 

-Gi-1 <^j, Bi 
. .-1 



[fc] 

-Buj 



0. 

[k] 



sequences {A} and {G} are computed by the formulae: 
(2.10) 



fAj+i = gj-PjAj Vj, A2 = gi, i=2,...,m, if A^t^ for all 2<i<m; 
I if Aj = for any i from (2<i <m),then Aj+i— is undefined, but Aj_|_2= ; 



(2.11) 



(Ik] Ik] Ik] Ik] 

Gi-i=?j -Tj+i Gi Pj+i,Gij^-i=gi^,i = ife-l,ifc-2,...,Zfe+i+l, if Giv^ 0; 

Ik] [k] 
if Gi=0 for any i from (Z^+i + l<i <Zfe — 1), then Gi-i— is undefined, 

Ik] 

hut Gi-^gi_i. 



[k] 



Ik] 



The structure elements (3^ and which determine the elements of submatrices B and 

[fc] 

their products H and H are computed [1] by the formulae: 



(2.12) 



' -p,A, \ if A, 7^ 0, [fc] 

-Pi: and pi+i = -Pi+iC^i, /3i+r 
. ifA, = 0; 



-r,+i Q , if G,7^ 0, 

[k] 

-ri+i, and /3i= -r^a;^, 

[fc] 

if G,= 0; 



(2.13) 



\(3i-pi-v --Pj+i, iij<i, 




[fc] [fcj [fc] 
1, if ^ > j; 



Proof. Let the system C-^X = Y, according to the theorem condition, be ill-posed 
but non-degenerate. Then its solution X+ with the properties given in the theorem does 
theoretically exist and it is unique. Let us show that it can numerically be obtained by 
the method (2.1) (2.13) called in [1] the critical-component method. To this end, we 
verify first that to the solution X+ there corresponds the following generahzed LDR [1] 
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decomposition of the matrix C3 (1.4): 

1 



(2.14) Cs^LDR- 



'2] [2] 
(Pl^+^l2l3+i)-(Pi^+^l2l2^ 1 



X 



X 



9i r2 
P2 Q2 ^3 



Na„+i+i=i) 



'?/3 + l^«3+2 
f'«3+25'i3+2^Z3+3 



PI2 %. 



[2] 



«3+l 



[1] 
C 



^2+1 
(Zl=m) 



P;2+2^'2+2^'2+3 



(Siirfin+i) 

[n] 



(Bi3+Uf(2+i) 



^^'2'J'2+l 
1 



where it is assumed that tridiagonal matrices c\'^~^^~^^ {k — 1,2, ...,n;li — ■m,ln+i + l — 1) 
are well-posed and their first diagonal elements are denoted by 



(2.14)' 



[k+i] 



Pik+1+1 B 1^+11^+1^1^, j_+i,k - 1,2,.., n 1, 



[k+i] 



where {Pi,Qi,ri}iJj^^^^+i are elements of the initial matrix C3 (1.4), B i^^^^j (j — 

[fc+i] 

Ik+i, Ik+i - 1, lk+2 + 1) are the last rows and B {i = k+i, k+i - 1, lk+2 + 1) are 

the last columns of matrices, inverse of the matrices ^ C 

[k] [k] 
with Bij (2.8) since they are elements of rectangular submatrices B (2.8). From the as- 



computed in accordance 



sumptions for C3 being nonsingular and for square matrices C/'^^^^^ 



being well-posed it 



follows that the LDR decomposition (2.14) is unique and stable to errors {h,ei,eo). 
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And the matrix B 



o 

— C3 can uniquely be represented in the form {B — {E + Q) B): 



(2.15) B^{B^ 



P2 I2 ^3 



P;3+2%+2'^i3+3 



PI2 %. 



,[2] 



'3+1^,-1 



Pm 



X 



X 



[3] 



("^(3+1^(3(4+1) 



[3] 



(-^,3 + 16,3,3) 



[2] [3] [2] [3] [2] [2] 

(P(2+iSi2'3+^'3+3^'3'4+i^"'^i2+i^'2'3+i^'3+i^'3'3^''"^'2+^'2'3+i^"'''"^'2+^'2'2^ 1 



+ (0 = 



[2] 



(-BiirTin+i) ^'"^ ('^"n''in+i'^'3+"2'''2+i'' 



U---U (,Sznirf(n+i^'3+"fi2+i^ 



[2] 



(~^(3+iifi2+i) 



, [2] 

'■~-^'2'i''2+l 

"0 



B 



) + 



Schematically, the matrix B can be represented as follows: 




li=m 



(2.15)' 



^3 + 1 



l\=m 



t 



[2] 

B 



[1] 
B 



-B . 



Representation (2.15) is easily established by a direct verification of the matrix equalities 

[k] 

C^B = E = BCs, with representations (2.8) ^ (2.13) taken into account for elements Bij 

o 

of the matrix B and decompositions of C3 and B given by (2.14) and (2.15). 

Now using representation (2.15) for B, we obtain components of the vector X+ in 
form (2.2) and (2.5). Prom (2.15) it follows that X+ can be written in the form 

(2.16) x+ = (E + n) x=x +nx, 

o 

where the vector X looks as follows 

(2.17) X=iB=C3~')Y 

and is a unique, stable to errors {h,S) and (£i,£o)) solution of the well-posed system of 
linear algebraic equations 



(2.18) (C3= 



Ql ^2 
P2 Q2 ^3 



W(/„+i+l=l) 
In 



9/3+l^«3+2 

PI2 %J 



[2] 

--C 



13+1 



c 



ill = rn) 







X- 



X, 



[2] 

m 

X 



12+1 



[1] 

x^ 



Vi 

Vln. 
1/(3+1 

yi2 , 
yi2+i 



which differs from the initial system CsX = y by the change of the corresponding off- 
diagonal elements to zeros and of diagonal elements q to elements q calculated by formulae 
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(2.14)'. Here the vector X includes components given by sums (2.2), which follows from 

o 

the representation B (2.15) and (2.8). 

For the matrix Q (2.15) we can write the following decomposition 



(2.19) Q = 







i-BurTin+i') 

[n] 









, [21 , 



. [21 , 





[2] 



, [2] 
'1 



Then for components of the vector il-X, denoted as vector 0, we obtain the exphcit form 

o [n] [n] [2] 

(2.20) [(- Sl/„r;„+i)x^+i, (- Blnlnfln+^^ln+l\ -A (" ^ i3+l'2'^i2+l)^S+l' - 



, [2] 



+ 



_ l;O,...,oF = 0. 

Consequently, components of the vector (f) are also calculated by formulae (2.2). 

As a result, wc have established that if is a normal (pseudo) solution of the system 
C3X — Y with the properties given by the theorem, then to it there correspond consistent 
with it decompositions (2.14) and (2.15) for matrix C3 and its (pseudo) inverse matrix 

B = C^. In this case, representations for X+ ^X +ftX (2.16) and +ft B (2.15) 

being consistent with each other (since they arc calculated with the same matrix Q (2.15)) 
are unique and stable to small errors h,6) and (£i,£o) in view of decompositions of C3 

o 

(2.14) and (2.15) being unique. Stabihty is a consequence of the matrix C3 (2.18) 
being well-posed. 

Now let us show that if the numerical solution of the system C3X — Y is obtained by 
the method*^ (2-1) -r- (2.13), then it is minimal in norm and provides a minimum of the 
discrepancy norm. Indeed, let X~^ is determined in the form (2.1) (2.13). Then from 
(2.5) it follows that the vector X~^ can be represented as a sum of two vectors. 



(2.21) X+ ^X +0, 

*)As we can see, this method includes the algorithm and criterion (2.3) (2.4) of separation of well- 
posed subspaccs and, respectively, the procedure of numerical finding of . It is to be kept in mind 
that the quantities (2.4) obey the inequalities |$^| < |Aj|, where A^- is a discrepancy. As a matter 

[k] [k] [k] 

of fact, \<^>j\ = \\yj\ - \yj - [y^ - (pj x j_i + q. x j + Vj+i x j+i)]|| = \\y.\ - \y. - A^.|| < |A^.|, since 
||y,|-|t/,-A,||<(|y,-(y, + A,)| = |A,|). 
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whose components are determined according to (2.2). The vectors X and 4> consist of n 
subvectors of proper dimensions, i.e., 

o_ [n] [n] [k] [k] [2] [2] [1] _ ^ [1] _ 

X— [[Xi, Xi^), ...,[X Xi^), ...,(^a;;g_|_i, Xi^), ...,[Xi^^i— X^^+i, Xm—x^\ , 

[n] [n] [k] [k] [2] [2] [1] [1] 

0=[(01, 0i„), i<Plk+i+l^ <^ife)' (<^i3+l' ^h)^ -' (<^i2+l^ 0' '^"^^ ' 

corresponding to well-posed n-subspaces that are separated in accordance with the crite- 
rion (2.3), (2.4). 

^From (2.1) ~ (2.5) we have that to the solution X+ there corresponds the decompo- 
sition of i? = of the form (2.15), which results in the representation for C3 (1.4) of 
form (2.14). Then, as mentioned above, the solution X~^ is written in form (2.16), where 

o 

X is in a unique way represented in form (2.17), (2.18). Owing to system (2.18) being 

o 

well-posed, which results from criterion (2.3) and (2.4), the vector X is unique, obeys the 

o 

condition min 1 1 X 1 1, and is stable to small errors {h, 5) and (^i, ^o)- From the uniqueness 
of matrix Q (2.15) that contains the last columns of matrices, inverse of the well-posed 

matrices C;^'*'^ , and from (2.16) and (2.15) it follows that X+ and {B = C^) are unique 
and minimal in norm. 

Let us now show that the vector X'^ determined by formulae (2.1) (2.5) satisfies the 
condition of minimum of the discrepancy norm (min HCaX"*" — y||). Taking advantage of 
the representation of X~^ (2.15), we get 

\\C3X+-Y\\ = WCsiE + Q) 1 = II CsX -Y\\, 

o 00-0^ o 

where X= {B=C^ )Y , , and B are, respectively, defined by (2.18) and (2.15). Owing 

00 - 

to the system (2.16) being well-posed, the minimum min || C3X —Y\ \ and, consequently, 
min ||C'3X+ —Y\ \ are attainable. So, the theorem is proved. 

Corollary. The norm of discrepancy HCaX^ — y|| complies with the following esti- 
mate: 
(2.22) 

[k] 

IC3X+- y||^ < siTp-f max \y-\ + A, where r = max (|p.|, |fj),p = max(| Bij |), 

l<i<m 2<i<m t,3',k 



1-\ T. hih - h+i), h = m, In+i = 0, < A < /i| |X+| I + 5. 

V '^=^ 

Proof. In view of all said above, we have 

||C'3X+-F|| = 11^3(^ + 1^) X -Y\\ = II C3X -Y\\. 

o o ~ o o ~ 

Since the system CiX= Y is well-posed, the Euclidean norm of errors || C3X — ^||e can 
be estimated by using the known results [9]: 

(2.23) \\C3X+-Y\\e ^WCsX -Y\\e < ^f{m)e,\\ C3 \\e\\ X We- 

However, in the case of the method considered above, this estimate turns out 
to be excessive. Actually, performing obvious transformations and making use 
of the definition of matrix norms consistent with the corresponding vector norms. 
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we get 
(2.24) 



CsX -Y\ 



o o - 



CsBY-YW = \\{C3B -E)Y\\ < \\ C3B -E\\ ■ \\Y\\. 



Next, we estimate the norm of matrix discrepancy || C3B —E\\, using the exphcit form 

00 o 

of C3 (2.18) and B (2.15), as well as the condition for the matrix C3 being well-posed. 

o o 

Taking account of the explicit form of elements of the matrix {CsB —E) and introducing 

o o 

the notation Vijiu — (uij) —C3B —E) for them, we can write the system of scalar identities 



(2.25) 



[fc] [k] [k] 

Pi Bi-ij +qi Bij +ri+i Bi+ij= i^ij, h+i + ^ < i < j < k, 

[k] [k] [k] , X , 

Pi Bi-ii +qi Bii +rj+i Bi+ii= 1 + ua, k+i + 1 < (« = j) < h, 
[k] [k] [k] 

I Pi Bi-ij +qi Bij +ri+i Bi+ij= i^ij, 1 < J < i < k- 



[k] 



Hereafter, k = 1,2, ...,n;li = m, = 0. Utilizing the representations for Bij (2.8) -j- 
(2.13), we write the system of identities (2.25) either in the form 



(2.26) 



+ ^i+i] Bi+ii+i n f3i= ^ij, k+i + 1 < « < j < 4, 

[k] [k] 

(Ai+i+ Gi-i -Qi) Bii= 1 + lyu, k+i + 1 < (« = j) < h, 

[k] [k] . ^ [k] i [k] 

[Gi-1 (- Gi\Pi) + Pi] Bi-ii-i n Pi= i^ij,l<J <i< h, 



or in the form 



(2.27) 



[k] 



(1 - Ai+iAi+i)ri+i Bi+lj= lk+1 + I <i< j < k, 

[k] ^ [k] 

(Ai+i+ Gi-i -Qi) Bii= 1 + Vii, k+i + 1 < (i = j) < Ik, 
[k] [k] [k] 
[ (1- Gi-iGi-i)Pi Bi-ij= i^ij, 1< j <i< Ik- 
Let us now estimate (2.27)i), (2.27)2) and (2.27)3); we have 

[k] 



(2.28) < 



[k] 



1(1- Ai+iAi^i)fi^^ Bi+ij= Uij\ < \1 - (1 ± ei)\ max \f-_^^ \ max | Bi+ij \ , 

!<«< m— 1 l<t<j<m—l 

[k] ^ [k] [k] 

|(Ai+i+ Gi-i -qi)Bii -1 = i/jjl < 1(1 ±£1) - 1| max \ Bii\, 



l<i<m 



[k] [k] [k] ^ [k] 

[I- Gi-iGi-i)PiBi-ij= i^ijl < |1 - (1 ±£i)| max |pj max | B 



2<i<h 



i-lj 



With estimates (2.28), we obtain \\C3X+ - Y\\^ = IKCaS -E)Y\\^ < \\ C3B 



-E\ 



M 



1^1 loo < W E h{lk - Wi) max|i/y|||y||oo < £iTp7 max |%|, where T,p,7 are de 



k=l 



l<i<m 



fined by (2.22). Here we took advantage of the condition of consistency of vector norms 
||C'3X+ — y||oo = niax \(CsX'^ — Y)i\ and ||l^||oo = max \yA with the M-norm of the 

l<i<m l<i<m 

. o o 

M = vm^ max \iC3B —E)ij\. The validity of in- 

l<i,j<m 

o o 

equality (2.22) is established. Since the Euchdean norm of the matrix || C3B —E\\e is 



matrix {C3B -E), i.e. || C3B -E\ 
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consistent only with the vector Euchdean norm | \Y\\e, instead of (2.22), one can, by anal- 

~ ^ o o — — 

ogous arguments, obtain the estimate IIC3X+ — yjj^; = IKCaS —E)Y\\e < eiTp^\\Y\\E, 
where T,p,'-f are defined by (2.22). 

Remark 1. To save the volume of publication, we do not present the method of 
solution of system (1.3) with the two-diagonal matrix C2 (1-4). It is expounded in detail 
in ref.[3] and it is shown there that it results from the method (2.1) (2.13). 

The estimate (2.22) for system (1.3) acquires the following form: 
(2.29) 

\\C2X+ -Y\ \ < II C2B -E\ \ • ||y|| < £iTp7 max |yj + A, where f = max(|fj), 

l<i<m 2<i<m 



[k] n 

p = max{\ Bij\),^= Jl/2 E {h - h+i), h = m, l^+i = 0, < A < /i| |X+| | + 5. 

[fc] 

Here Bij are elements of upper triangular matrices, inverse of well-posed two-diagonal 

rt''] n'fc-H+1 

matrices [6*2]/^ • 

Remark 2. Note that due to orthogonality of matrices P and Q in transformations 
(1.5), the following estimates take place for the norms of discrepancy \\AZ^ ~ -^11- 



(2.30) 



I |AZ+ - -F| I < £ifp7 max |7/.| + A, \i A = A^, 

l<i<m 

I \AZ+ -F\\< eiTpi max |y^| + A, li A ^ A^, 

l<i<m 



where f, p, 7 and f, p, 7 are defined in analogy with (2.22) and (2.29), < A < h\ \Z'^\ \ +6. 

Remark 3. The above estimates (2.22), (2.29) and (2.30) can also be used for prob- 
lems of inversion i.e., C^C^ = E, C2C2 = E, AA^ = E: the matrices C^, C2 and A'^ are 
to be obtained by solving the matrix system of equations 

C3C3+ = E, C2C+ = E and AA+ = E. 

by the critical-component method. In the case when systems (1-1), (1-2) and (1.3) are 
ill-posed, one should not take, as C^, C2 and the corresponding matrices obtained by 
the critical- component method in solving these systems of equations with a given right- 
hand side. The reason is that the norms of matrices C3", C2 are consistent with the 

norms of concrete vectors X^, X^, Y, Y and F. 

Remark 4. The theorem is formulated under the assumption detC* 7^ 0. Let us re- 
move this restriction. The critical-component method does not explicitly use the quantity 
detC. Rather, it is based on the processes (2.10) and (2.11) for computing elements of 
m-dimensional vectors {A, G}. As established in ref. [14], if detC = 0, then components 
of these vectors get into one of the following three situations: either A^+i = and Gq = 0, 
or Aj = and Gj = 0, or [(Aj = and AjAj+i = 0) or (G., = and GiGi_i = 0)]. In 
this case we replace some zero quantities by the quantity o{6i). This does not essen- 
tially impair the quality of solution, since such perturbations can already be present in 
these quantities. Gonsequently, one may consider the critical-component method to be 
applicable for any value of det (7, including det C — 0. 

3. Results of numerical; experiments and their analysis 

In this section, we discuss the results of numerical experiments performed in the com- 
puter arithmetic with double accuracy (£1 = 2~^^ Ki 2,2 • 10~^^) for computing basic 
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numerical characteristics of the solutions X of systems WX — W : Ci\C'i\A ^ 
A^\A = . Let us first explain the notation and abbreviations adopted in Tables 
1 -i- 12 : 5^-* — the relative error of X^"*) — the obtained numerical solution of 
system W^M^^^") = F^™) {W^""^ : ]Ct^ ; A^"^^ + (^M^t.^M _ (^M)T _ 

the above indicated types of matrices, X^"*) — the exact solution, m — the order of 
the system under consideration); //(VF*^"*)) = cond(VF^"*^) — the condition number of 
iy(m). ) ^ com.time(sec.) — the time of computing solutions X^'"), ^['"^ 4""^ 



— the lower and upper bounds t^^-* i.e.*) 
where HX^^^H^; are norms of approximate and exact solutions; Sl - 

■'1=1 



N, , , A^, , , N, N. 

1 



^1=1 U=i ^ U=i ^1=1 



Nj Nj 

jjL{W) — . S (Ai(VF("*)))/, i(sec.) = jr. J2it^"^'')i are arithmetic means of the character- 

^1=1 U=i 

_ s 

istics listed above, Ni = J2 ^j: where A^^ is the number of examples of a given type, 

i=i 

s - the number of examples in a Table, i - the number of a Table; MCS and MCC — 
our programs DCSOL (access through www http : //cv . j inr . ru/lcta/ sap/ lib/f 499 . f ) 
from library LIBJINR [17] (algorithms of the critical-component method [1 -j- 3]); GS 
— programs DBEQN and DEQN from library CERNLIB [21] (a modified algorithm of the 
Gauss exclusion method); OSM — program DTSYS from library LIBJINR [19] (algorithm 
of nonmonotone orthogonal run); QR — programs FOIAXF from library NAGLIB [20] (al- 
gorithms of QR — method); SVD — subprogram- function PSOL from library LINA [7] 
(algorithm of the singular-expansion method with the use of exhaustion); TRM — sub- 
program SLAY from library LIBJINR [19] (algorithms of the Tikhonov regularization 

10 _ 

method). We have solved N ^ 278 {N ^ E A^j, i 4 and i 8) different sys- 

i=l 

tems of linear algebraic equations at different orders**^ {k = 1, 2, 3, 4, 4', 5), which 
is presented below. A system of the type C2X — Y. Examples 1 -j- 5 from [2,3,16] 
(see §5 Appendix), {mi : 10; 20. m2 : 3. : 5; 10; 15. : 10; 20. Here 
1 < ^(C(™fc)) < 1/^}, {nil ■■ 30; 40; 50. ms : 4; 5; 6. mg : 20; 25; 30; 35. : 30; 40. 
Here 1/^ < /i(C(™fe)) < 1/si}, {mi : 60; 70; ...; 100; 150; 200. : 8; 9; 10. rus : 
40; 45. rui : 5; 6;...; 18. : 50; 60; ...; 100; 150; 200; 300; 400; 500. Here 1/ei < fx{C^"'k))}. 
A system of the type C3X — Y. Examples''^*) . 6 10 from [2,3,16] (see §5 Appendix). 



*^The left-hand side of inequality (3.1) is a property [10] of the norm || • ||, and the right-hand side 
is obtained by using the exact solution X = W~^Y and equality W~^W = E. We have {5m = \\X — 
= \\X - W-^Y\\/\\X\\ = \\W-^{WX - < {\\W-^\\ ■ \\WX - Y\\/\\X\\ = 6r). Note 

that in practice, inequalities (3.1) can be broken (see, for instance. Table 7). This occurs when calculating 
WX — Y. In this case, the solution X can be surely considered acceptable. 

lower index k of order indicates the number of an example from the set of given- type examples. 
■f*'Example 4' is example 4 from [2] (system 9, see §5 Appendix), but with =0,00000001. 
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{mi : 10; 20; ...; 100; 150; 200; 300; ...; 900. rui : 10; 20; ...; 100; 150; 200; 300; ...; 900 

mg : 10; 21; 3 0; 40; 51; 60; 70; 81; 90; 100; 151; 201; 300; 400; 501; 600; 700; 801; 900 
rui : 10; 30; 40; 60; 70; 90; 100; 300; 400; 600; 700; 900. m^, : 10; 30; 40; 60; 70; 90; 100 
300; 400; 600; 700; 900. mg : 10; 20; ...; 100. Here 1 < /x(C(™fe)) < 1/^}, {m4 : 20; 50 
80; 150; 200; 500; 800. mg : 150; 200. Here 1/^ < //(C("^fe)) < 1/si}, {m^' : 20; 50 
80; 150; 200; 500; 800. Here 1/si < /x(C('"fe))}. A system of the type AX^Y {A^ A^) 
Examples 11 ^ 15 from [2,3,16] (see §5 Appendix), {mi : 10; 20; ...; 100; 150; 200; 300 
400; 500. m2 : 5. mg : 10; 20; ...; 100; 150; 200; 250; 300; 400; 500. : 3; 4; ...; 10 
ms : 5; 6. Here 1 < /i(C(™fc)) < 1/^}, {ms : 6; 7; ...; 10. : 11; 12; ...; 17. mg 
7; 8; ...;11. Here 1/^ < /x(C('"fe)) < l/^i}, {m-i : 11; 12. : 18. mg : 12. Here 
l/£i < fiiC^'^k))}. A system of the type AX ^ Y {A ^ A^). Examples 16 ^ 20 from 
[2,3,16] (see §5 Appendix), {mi : 10; 20; ...; 100; 150; 200; 300. mg : 5; 6. mg : 5; 6. m4 : 
10; 20; 100; 150; 200; 300; 400; 500. Here 1 < ^(C('"fc)) < 1/^}, {ms : 7; 8; ...; 11. 
ma : 7; 8; ...; 11. Here 1/^ < /x(C("*fc)) < l/^i}, {ms : 12. ma : 12. mg : 5; 10; 15; 20; 
25; 30; 35. Here 1/ei < fx{C^'^k))}. 

Below, in Tables 1 3, 5 7 and 9 11, we report the obtained numerical values of 
the indicated characteristics S^, Sfi, 5^, Sx, of approximate X and exact X solutions of 
systems WX = Y {W : C2; C3; A^ A^;A = A^) at 1 < i^{W) < 1/^ — being well- 
posed, 1/v^ < IJ'{W) < 1/ei — being ill-posed, 1/ei < IJ,{W) — being pathologically 
ill-posed of these systems, respectively. In Tables 4, 8 and 12, we present averaged results 
of Tables 1 3, 5 7 and 9 11. Note also that when l/si < //(VF), the subprogram 
SVD stops to work producing information INF= —1. The program TRM docs not work [19], 
when m > 100. Tables 9, 12 do not contain (*****) values of t{.) above 100 sec. 

For an easier apprehension of the calculation results reported in Tables 4, 8 and 12, 
we plot their "graphic images-figures" . 

Remcirk 5. In Tables 4, 8 and 12 (as well as in Figures 1 9) , we also present 



the averaged results of computations by subprograms MCS when W = {A = A^) and 



DSM I when W = C3, which is to be kept in mind when analyzing the above Tables and 
Figures. Explanations to some Tables and Figures: on the horizontal axis of Figs. 1 -j- 9, 
in units 10* where values of order t are indicated below the axis, approximate values of 
quantities 6l{W), Sm{W) and Sr{W) from Tables 4,8 and 12 are plotted. On the vertical 
axis of these Figures, we point out the names of programs through which those values 
have been obtained. The horizontal axis of Figs. 11 ^ 13 represents relative errors of the 
r.h.s. of a system with the Hilbert matrix (see §5 of Appendix, example 17), whereas on 
their vertical axes, we plot the corresponding relative errors found by various programs. 
For names of programs, see the notation. Table 13 contains the numerical results that 
arc drawn in Figs. 11 -^ 13. Note that: < 6y>=< n^Trr > ^ the mean value of the 



relative error of perturbation of the r.h.s. of the system; < Sx>=< ' ^ > ^ exact mean 
values corresponding to < 5y > . Numbers in Table 13 written in line with a program are 
average values really obtained by this program for <Sx> . In Fig. 10, the matrix, inverse 
of the Hilbert matrix, of order m = 14 is graphically shown. Along the axis Z, the values 
of elements of this inverse matrix are indicated. As a result, its complicated structure is 
easily visualized. Also, numbers of subspaces separated by the critical-component method 
are given in the Figure . 
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Table 1 (0.173E02< fi(W) < 0.547E08, 3< m <100, fl{W) = 0.533E07, iVi = 117) 





PR. 


t(scc) 


5l 




Sr 




Sx 




MCC 


0.0004 


0.238E-11 


0.301E-11 


0.470E-10 


0.136E01 


0.136E01 


C2X = Y 


GS 


0.0005 


0.238E-11 


0.301E-11 


0.119E-08 


0.136E01 


0.136E01 


Ni = 8 


SVD 


0.1568 


0.261E-11 


0.335E-11 


0.723E-08 


0.136E01 


0.136E01 




QR 


0.0036 


0.113E-10 


0.175E-10 


0.405E-08 


0.136E01 


0.136E01 


0.576E07 


TRM 


0.0366 


0.448E-05 


0.546E-05 


0.575E-05 


0.136E01 


0.136E01 




MCC 


0.0051 


0.839E-13 


0.868E-13 


0.137E-12 


0.609E01 


0.609E01 


C3X = Y 


GS 


0.0028 


0.994E-13 


0.153E-11 


0.289E-10 


0.609E01 


0.609E01 


N2 = 54 


QR 


6.9793 


0.187E-11 


0.302E-10 


0.294E-08 


0.609E01 


0.609E01 




SVD 


0.3091 


0.609E-12 


0.505E-10 


0.704E-09 


0.609E01 


0.609E01 


0.106^07 


TRM 


5.3795 


0.238E-07 


0.212E-05 


0.212E-05 


0.609E01 


0.609E01 




OSM 


0.0039 


0.561E-01 


0.164E00 


0.850E04 


0.655E01 


0.609E01 


AX = Y 


QR 


0.1776 


0.526E-11 


O.lOlE-09 


0.122E-07 


0.254E01 


0.254E01 


A^A^ 


MCC 


1.0925 


0.973E-11 


0.117E-09 


0.461E-08 


0.254E01 


0.254E01 


N3 = 31 


SVD 


2.8117 


0.973E-11 


0.117E-09 


0.683E-08 


0.254E01 


0.254E01 


Hi A) = 


GS 


0.1091 


0.752E-11 


0.144E-09 


0.256E-08 


0.254E01 


0.254E01 


0.549E07 


TRM 


2.9615 


0.364E-04 


0.618E-04 


0.639E-04 


0.254E01 


0.254E01 




MCC 


0.9636 


0.163E-10 


0.951E-10 


0.178E-07 


0.434E01 


0.434E01 


AX = Y 


SVD 


2.5378 


0.163E-10 


0.951E-10 


0.307E-07 


0.434E01 


0.434E01 


A = A^ 


QR 


0.1564 


0.129E-10 


0.160E-09 


0.383E-07 


0.434E01 


0.434E01 


N4 = 24 


GS 


0.0977 


0.373E-10 


0.179E-09 


0.230E-07 


0.434E01 


0.434E01 


fl{A) = 


MCS 


0.9367 


0.153E-09 


0.429E-09 


0.171E-07 


0.434E01 


0.434E01 


0.900E07 


TRM 


2.5976 


0.603E-08 


0.328E-07 


0.453E-07 


0.434E01 


0.434E01 


Table 2 (0.968E08< fi{W) <0.399E16, 6< m <^ 


^0, li{W) = 0.365E15, N2 = 43) 




PR. 


t(sec) 


Sl 




6r 


6*- 


6x 




MCC 


0.0008 


0.191E-04 


0.353E-04 


0.289E-02 


0.180E01 


0.180E01 


C2X = Y 


GS 


0.0009 


0.191E-04 


0.353E-04 


0.517E-02 


0.180E01 


0.180E01 


Ni = 13 


SVD 


0.7276 


0.598E-04 


O.llOE-03 


0.259E00 


0.180E01 


0.180E01 


r'(C2) = 


QR 


0.0200 


0.954E-03 


0.221E-02 


0.276E00 


0.180E01 


0.180E01 


0.351E15 


TRM 


0.3108 


0.208E12 


0.208E12 


0.210E12 


0.131E13 


0.180E01 




MCC 


0.0067 


0.638E-08 


0.240E-04 


0.277E-02 


0.754E01 


0.754E01 


C3X = Y 


GS 


0.0035 


0.672E-06 


0.213E-03 


0.313E-02 


0.754E01 


0.754E01 


N2 = 3 


SVD 


8.6764 


0.335E-05 


0.463E-03 


0.453E-01 


0.754E01 


0.754E01 


fiiCs) = 


QR 


0.4031 


0.102E-04 


0.962E-03 


0.960E-02 


0.754E01 


0.754E01 


0.586E15 


DSM 


0.0048 


0.468E-01 


0.141E00 


0.136E11 


0.793E01 


0.754E01 




TRM 


6.6831 


0.148E14 


0.148E14 


0.562E14 


0.125E15 


0.754E01 


AX=Y 


MCC 


1.0966 


0.200E-03 


0.418E-03 


0.239E00 


0.254E01 


0.254E01 


A^A^ 


SVD 


2.8427 


0.200E-03 


0.418E-03 


0.255E00 


0.254E01 


0.254E01 


iV3 = 17 


GS 


0.1095 


0.278E-03 


0.514E-03 


0.473E00 


0.255E01 


0.254E01 


fl{A) = 


QR 


0.1784 


0.338E-03 


0.615E-03 


0.316E00 


0.255E01 


0.254E01 


0.271E15 


TRM 


2.9711 


0.133E12 


0.133E12 


0.261E12 


0.167E12 


0.254E01 




MCC 


0.9660 


0.654E-06 


0.161E-04 


0.215E00 


0.435E01 


0.435E01 


AX = Y 


SVD 


2.5544 


0.654E-06 


0.161E-04 


0.467E00 


0.435E01 


0.435E01 


A = A^ 


MCS 


0.9389 


0.119E-05 


0.548E-04 


0.189E00 


0.435E01 


0.435E01 


Ni = 10 


GS 


0.0979 


0.166E-05 


0.828E-04 


0.452E00 


0.435E01 


0.435E01 


p.{A) = 


QR 


0.1568 


0.109E-05 


0.976E-04 


0.734E00 


0.435E01 


0.435E01 


0.253E15 


TRM 


2.6032 


0.435E05 


0.435E05 


0.460E05 


0.543E05 


0.435E01 
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Tabic 3 {fi{W) >0.450E16, 5< m <80, p,(W) >0.450E16, N3 = 46) 





rix. 


t(scc) 








^\ 






MCC 


0.0010 


0.133E00 


0.179E00 


0.952E02 


0.259E01 


0.253E01 


C2X = Y 


GS 


0.0012 


0.482E00 


0.103E01 


0.924E03 


0.378E01 


0.253E01 


Ni = 30 


QR 


r\ r\c> ''7r\ 

0.0879 


0.274E02 


0.276E02 


0.211E07 


0.182E03 


0.253E01 


fi[C2) > 


iHM 


i.byoo 


U.bOziliOO 


U.dOzIIjOO 


0.145E75 


0.602E56 


0.253E01 


0.450i!jlo 


SVD 






INF 


= - 1 








MCC 


0.0069 


0.298E00 


0.828E00 


0.126E02 


0.887E01 


0.683E01 


C3X = Y 


QR 


0.1659 


0.724E00 


0.127E01 


0.398E02 


0.129E02 


0.683E01 


iV2 = 3 


GS 


0.0022 


0.174E02 


0.184E02 


0.940E03 


0.142E03 


0.683E01 


/^(<-^3) > 


□SM 


U.U(Jo4 






0.140E18 


0.178E03 


0.683E01 


0.450r!jlD 


TRM 


3.1464 


0.777E15 


0.777E15 


0.611E16 


0.430E16 


0.683E01 




SVD 






INF 


= - 1 






AX = Y 


MCC 


0.0179 


0.109E-02 


0.333E-01 


0.106E03 


0.856E00 


0.856E00 


A 1 aT 


GS 


0.0022 


0.330E-02 


0.681E-01 


0.193E03 


0.861E00 


0.856E00 


N3 = 4 


QR 


0.0040 


0.610E-02 


0.895E-01 


0.209E03 


0.861E00 


0.856E00 


m(^) > 


IRM 


U.U4zl 


U.ibDilii4 


U.iDDilii4 


0.278E14 


0.208E14 


0.856E00 


0.450E16 


SVD 






INF 


= - 1 








MCC 


0.0479 


0.428E-01 


0.365E00 


0.139E03 


0.994E00 


0.989E00 


AX = Y 


MCS 


0.0439 


0.197E00 


0.545E00 


0.149E03 


0.107E01 


0.989E00 


A = A^ 


GS 


0.0051 


0.104E05 


0.104E05 


0.297E22 


0.483E04 


0.989E00 


iV4 = 9 


QR 


0.0090 


0.261E12 


0.261E12 


0.332E25 


0.123E12 


0.989E00 


n{A)> 


TRM 


0.1165 


0.457E16 


0.457E16 


0.593E18 


0.213E16 


0.989E00 


0.450E16 


SVD 






INF 


= - 1 






Table 4 (Mean values of characteristics of Tables 1 




, iV4 = 206) 




PR. 


t(sec) 


Sl 




Sr 




Sx 




MCC 


0.5154 


0.712E-11 


0.538E-10 


0.561E-08 


0.358E01 


0.358E01 


m) = 


SVD 


3.1214 


0.763E-11 


0.614E-10 


0.119E-07 


0.358E01 


0.358E01 


0.533E07 


GS 


0.0525 


0.118E-10 


0.819E-10 


0.669E-08 


0.358E01 


0.358E01 


Ni = 117 


QR 


0.1617 


0.752E-11 


0.822E-10 


0.138E-07 


0.358E01 


0.358E01 






MCS 




0.9367 


0.153E-09 


0.429E-09 


0.171E-07 


0.434E01 


0.434E01 




TRM 


2.7438 


0.102E-04 


0.174E-04 


0.180E-04 


0.358E01 


0.358E01 






OSM 




0.0039 


0.561E-01 


0.164E00 


0.850E04 


0.655E01 


0.609E01 






MCS 




0.9389 


0.119E-05 


0.548E-04 


0.189E00 


0.435E01 


0.435E01 




MCC 


0.5175 


0.549E-04 


0.123E-03 


0.115E00 


0.406E01 


0.406E01 


m) = 


GS 


0.0529 


0.749E-04 


0.211E-03 


0.233E00 


0.406E01 


0.406E01 


0.365E15 


SVD 


3.7003 


0.660E-04 


0.252E-03 


0.257E00 


0.406E01 


0.406E01 


N2 = 43 


QR 


0.1896 


0.326E-03 


0.971E-03 


0.334E00 


0.406E01 


0.406E01 






OSM 




0.0048 


0.468E-01 


0.141E00 


0.136E11 


0.793E01 


0.754E01 




TRM 


3.1420 


0.379E13 


0.379E13 


0.142E14 


0.316E14 


0.406E01 




MCC 


0.0184 


0.119E00 


0.351E00 


0.882E02 


0.333E01 


0.280E01 


m)> 




MCS 




0.0439 


0.197E00 


0.545E00 


0.149E03 


0.107E01 


0.989E00 


0.450E16 


OSM 


0.0034 


0.221E02 


0.231E02 


0.140E18 


0.178E03 


0.683E01 


iV3 = 46 


GS 


0.0027 


0.260E04 


0.260E04 


0.742E21 


0.124E04 


0.280E01 




QR 


0.0667 


0.653E11 


0.653E11 


0.830E24 


0.308E11 


0.280E01 




TRM 


1.2501 


0.151E55 


0.151E55 


0.362E74 


0.151E56 


0.280E01 




SVD 






INF 


= - 1 
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Table 5 (0.250E03<;u(VF) <0.323E07, 150<m<200,Af(Pr)= 0.827E08, iV5= 16) 





PR. 


t(sec) 


5l 


Sm 


Sr 


Sl 


Sx 




MCC 


0.0133 


0.344E-12 


0.362E-12 


0.724E-12 


0.788E01 


0.788E01 


C3X = Y 


GS 


0.0065 


0.345E-12 


0.363E-12 


0.727E-12 


0.788E01 


0.788E01 


N2 = 8 


QR 


6.0108 


0.622E-12 


0.777E-12 


0.170E-10 


0.788E01 


0.788E01 


= 


SVD 


31.137 


0.138E-10 


0.104E-10 


0.274E-10 


0.788E01 


0.788E01 


0.202E05 


OSM 


0.0093 


0.306E00 


0.467E00 


0.114E01 


0.116E02 


0.788E01 


AX = Y 


MCC 


22.606 


0.331E-13 


0.311E-09 


0.120E-07 


0.719E01 


0.719E01 


A^A^ 


GS 


2.2529 


0.134E-13 


0.374E-09 


0.275E-07 


0.719E01 


0.719E01 


A^3 = 4 


SVD 


37.433 


0.331E-13 


0.436E-09 


0.189E-07 


0.719E01 


0.719E01 


jl{A) = 


QR 


3.6669 


0.236E-13 


0.617E-09 


0.306E-07 


0.719E01 


0.719E01 


0.123E07 
















AX = Y 


MCC 


10.312 


0.122E-11 


0.155E-09 


0.121E-06 


0.134E02 


0.134E02 


A = A^ 


SVD 


18.469 


0.122E-11 


0.155E-09 


0.364E-06 


0.134E02 


0.134E02 


Af4 = 4 


MCS 


10.181 


0.115E-11 


0.215E-09 


O.lOOE-06 


0.134E02 


0.134E02 


p.{A) = 


QR 


1.6768 


0.808E-12 


0.263E-09 


0.371E-06 


0.134E02 


0.134E02 


0.123E07 


GS 


1.0537 


0.218E-11 


0.336E-09 


0.230E-06 


0.134E02 


0.134E02 


Table 6 (0.657E1K /i((7 


3) <0.594E15, 150< m < 


200, mCs) -- 


= 0.308E15, 


^6 = 3) 




PR. 


t(sec) 


Sl 


Sm 


Sr 




Sx 




MCC 


0.0591 


0.117E-05 


0.133E-02 


0.207E-01 


0.136E02 


0.136E02 


C^X = Y 


SVD 


43.881 


0.488E-05 


0.179E-02 


0.171E01 


0.136E02 


0.136E02 


N2 = 3 


GS 


0.0075 


0.766E-05 


0.277E-02 


0.324E-01 


0.136E02 


0.136E02 


HC3) = 


QR 


4.5427 


0.114E-04 


0.336E-02 


0.229E00 


0.136E02 


0.136E02 


0.308E15 


OSM 


0.0109 


0.135E-02 


0.414E-01 


0.101E12 


0.138E02 


0.136E02 


Ta 


ble 7 {n{W) >0.450E16, 150< m <200, il{W) >0.450E16, Nj = 


5) 




PR. 


t(sec) 


Sl 


Sm 


Sr 


S*- 


Sx 


C2X = Y 


MCC 


0.0049 


0.615E-04 


0.113E-03 


O.OOOEOO 


0.719E01 


0.719E01 


Ni=4 


GS 


0.0060 


0.615E-04 


0.113E-03 


O.OOOEOO 


0.719E01 


0.719E01 


-KC2) > 


QR 


3.3505 


0.401E02 


0.408E02 


0.128E06 


0.653E02 


0.719E01 


0.450E16 


SVD 






INF 


= - 1 








MCC 


0.0278 


0.292E+00 


0.819E+00 


0.258E02 


0.183E02 


0.141E02 


QiX = Y 


QR 


5.3062 


0.551E+00 


0.837E+00 


0.446E02 


0.189E02 


0.141E02 


N2 = 1 


GS 


0.0079 


0.719E+02 


0.728E+02 


0.649E04 


0.103E04 


0.141E02 




OSM 


0.0133 


0.882E+02 


0.892E+02 


0.571E18 


0.126E04 


0.141E02 


0.450E16 


SVD 






INF 


= - 1 






Table 8 


(Mean values of characteristics of Tables 5 


7, iVg = 24) 




PR. 


t(sec) 


Sl 


Sm 


Sr 




Sx 




MCC 


10.977 


0.532E-12 


0.155E-09 


0.443E-07 


0.949E01 


0.949E01 


m) = 


SVD 


29.013 


0.502E-11 


0.200E-09 


0.128E-06 


0.949E01 


0.949E01 


0.827E06 


MCS 


10.181 


0.115E-11 


0.215E-09 


O.lOOE-06 


0.134E02 


0.134E02 


iVs = 16 


GS 


1.1044 


0.846E-12 


0.237E-09 


0.848E-07 


0.949E01 


0.949E01 




QR 


3.7848 


0.498E-12 


0.294E-09 


0.134E-06 


0.949E01 


0.949E01 




OSM 


0.0093 


0.306E00 


0.467E00 


0.114E01 


0.116E02 


0.788E01 




MCC 


0.0591 


0.175E-05 


0.133E-02 


0.207E-01 


0.136E02 


0.136E02 


^(Cs) = 


SVD 


43.881 


0.488E-05 


0.179E-02 


0.171E01 


0.136E02 


0.136E02 


0.308E15 


GS 


0.0075 


0.766E-05 


0.277E-02 


0.324E-01 


0.136E02 


0.136E02 


iVe = 3 


QR 


4.5427 


0.114E-04 


0.336E-02 


0.229E00 


0.136E02 


0.136E02 




OSM 


0.0109 


0.135E-02 


0.414E-01 


0.101E12 


0.138E02 


0.136E02 
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0.450E16 






MCC 
QR 
GS 




0.0164 
4.3284 
0.0069 


0.146E00 
0.203E02 
0.360E02 


0.410E00 
0.208E02 
0.364E02 


0.129E02 
0.640E05 
0.325E04 


0.127E02 
0.421E02 
0.519E03 


0.106E02 
0.106E02 
0.106E02 


iV7 = 5 






OSM 




0.0133 


0.882E02 


0.892E02 


0.571E18 


0.126E04 


0.141E02 








SVD 








INF 


= - 1 






Table 9 (0.498E03< ii{W) <0.318E08, 250< m <900, li{W) = 


= 0.675E07, 


Ng = 39) 




PR. 


t(sec) 


5r 


6m 


Sr 

11, 


"x 


5x 


C^X = Y 
HCs) = 
0.255E06 

iV2 = 31 


MCC 
GS 
OSM 


0.0524 
0.0230 
0.0324 


0.762E-14 
0.718E-14 
0.179E02 


0.403E-13 
0.444E-13 
0.186E02 


0.168E-11 
0.179E-11 
0.682E75 


0.999E01 
0.999E01 
0.618E02 


0.999E01 
0.999E01 
0.999E01 


AX = Y 
A + A^ 

A^3 = 4 




MCC 
GS 






12.156 
1.2256 


0.703E-12 
0.227E-11 


0.611E-08 
0.971E-08 


0.115E-05 
0.179E-05 


0.101E02 
0.101E02 


0.101E02 
0.101E02 


Af4 = 4 
A = A^ 

m) = 

0.999E07 




MCC 
MCS 
GS 








0.256E-11 
0.216E-12 

0.703E-12 


0.522E-09 
0.227E-08 

0.240E-08 


0.182E-05 
0.174E-05 

0.197E-05 


0.188E02 
0.188E02 

0.188E02 


0.188E02 
0.188E02 
0.188E02 


Table 10 (0.540E15< /i(C3) <0.637E15, m : 500,800, lL(C-i) = 


=0.589E15, iVio = 2) 




PR. 


t(sec) 


Sl 








5x 


C3X = Y 

HCs) = 
0.589E15 
N2 = 2 


MCC 
GS 
OSM 


0.0901 
0.0290 
0.0394 


0.161E-05 
0.294E-04 
0.150E-02 


0.174E-02 
0.765E-02 
0.915E-01 


0.415E-01 
0.666E-01 
0.290E75 


0.253E02 
0.253E02 
0.253E02 


0.253E02 
0.253E02 
0.253E02 


Table 11 {f-i{W) >0.450E16, 250< m < 800, 


Ji{W) > 0.450E16, A^ii 


= 7) 




PR. 


t(scc) 


Sl 










C2X = Y 

A(C2) > 
0.450E16 

iVi = 4 


MCC 
GS 


0.0088 
0.0115 


0.615E-04 
0.615E-04 


0.113E-03 
0.113E-03 


O.OOOEOO 
O.OOOEOO 


0.106E02 
0.106E02 


0.106E02 
0.106E02 


C3X = Y 
A(C3) > 

0.450E16 
iV2 = 3 


MCC 
GS 

OSM 


0.0583 
0.0192 

0.0300 


0.145E00 
0.167E03 
0.215E03 


0.409E00 
0.168E03 
0.216E03 


0.114E03 
0.190E06 
0.209E75 


0.250E02 
0.394E04 
0.500E04 


0.213E02 
0.213E02 
0.213E02 


Table 12 (Mean values of characteristics of Tables 9 -r- 


11, Ni2 = 48) 




PR. 


t(sec) 


6l 




5r 




5x 


fl{W) = 




MCC 








0.109E-11 


0.221E-08 


0.990E-06 


0.130E02 


0.130E02 


0.675E07 




MCS 








0.216E-12 


0.227E-08 


0.174E-05 


0.188E02 


0.188E02 


Ng = 39 






GS 








0.993E-12 


0.404E-08 


0.125E-05 


0.130E02 


0.130E02 






OSM 








0.179E02 


0.186E02 


0.682E75 


0.618E02 


0.999E01 


A(C3) = 
0.589E15 




MCC 
GS 






0.0901 
0.0290 


0.161E-05 
0.294E-04 


0.174E-02 
0.765E-02 


0.415E-01 
0.666E-01 


0.253E02 
0.253E02 


0.253E02 
0.253E02 


iVio = 2 




OSM 






0.0394 


0.150E-02 


0.915E-01 


0.290E75 


0.253E02 


0.253E02 


HW)> 

0.450E16 




MCC 
GS 






0.0335 
0.0153 


0.725E-01 
0.835E02 


0.205E00 
0.840E02 


0.570E02 
0.950E05 


0.178E02 
0.198E04 


0.159E02 
0.159E02 


iVii = 7 




OSM 






0.0300 


0.215E03 


0.216E03 


0.209E75 


0.500E04 


0.213E02 
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Table 13 





3 < m < 6, 




7< m< 11, 




12 < m < 13, 




0.524E03</x(A)<0.150E08 


0.475E09</i(A) <).518E15 




H{A) >0.450E16 




<5x>= 0.00, <Sy>= 0.00 


<Sx>=Om, <(5r>=0.00 


<5x>=om, <5y>= 


- n nn 


MCS 


0.384886772512E-10 


MCC 


0.64196E-04 


MCS 


0.182E01 




QR 


0.394340547190E-10 


SVD 


0.64196E-04 


MCC 


0.214E01 




GS 


0.432343309209E-10 


MCS 


0.23960E-03 


GS 


0.270E01 




MCC 


0.502178164479E-10 


GS 


0.23988E-03 


QR 


0.403E01 




SVD 


0.502178164479E-10 


QR 


0.44822E-03 


SVD 


INF= -1 




TRM 


0.334842097124E-07 


TRM 


0.58231E05 


TRM 


0.397E08 




<^^>=0.10, <(5y>= 0.054 


<,5x>=0.10, <(5r>= 0.047 


<,5x>=0.10, <dY>= 


n 04^ 


MCS 


0.999999999782E-01 


MCC 


0.99945E-01 


MCS 


0.125E01 




MCC 


0.999999999809E-01 


SVD 


0.99945E-01 


MCC 


0.277E01 




SVD 


0.999999999809E-01 


QR 


0.99982E-01 


GS 


0.333E01 




GS 


0.999999999948E-01 


MCS 


0.10005E00 


QR 


0.465E01 




QR 


0.999999999991E-01 


GS 


0.10005E00 


SVD 


INF= -1 




TRM 


0.100000344252E00 


TRM 


0.23802E06 


TRM 


0.442E10 




<Sx>=0.20, <(5y>= 0.107 


<Sx>=0-20, <(5y>= 0.094 


<,5x>=0.20, <5y>= 


riQn 


TRM 


0.199999724752E00 


QR 


0.20003E00 


MCC 


0.219E01 




MCC 


0.199999999988E00 


GS 


0.20004E00 


MCS 


0.323E01 




SVD 


0.199999999988E00 


MCS 


0.20009E00 


QR 


0.342E01 




MCS 


0.199999999994E00 


MCC 


0.20016E00 


GS 


0.352E01 




GS 


0.199999999997E00 


SVD 


0.20016E00 


SVD 


INF= -1 




QR 


0.199999999999E00 


TRM 


0.10208E06 


TRM 


0.139E11 




<,5x>=0.30, <5y>= 0.161 


<<^^>=0.30, <5y>= 0.142 


<,5x>=0.30, <5y>= 




MCC 


0.299999999983E00 


MCS 


0.29990E00 


MCS 


0.172E01 




SVD 


0.299999999983E00 


GS 


0.29991E00 


MCC 


0.279E01 




MCS 


0.299999999985E00 


QR 


0.29995E00 


GS 


0.350E01 




GS 


0.299999999992E00 


MCC 


0.29996E00 


QR 


0.465E01 




QR 


0.299999999996E00 


SVD 


0.29996E00 


SVD 


INF= -1 




TRM 


0.300000398356E00 


TRM 


0.60699E04 


TRM 


0.487E10 




<,5x>=0.39, <(5y>= 0.209 


<,5x>=0.39, <(5y>= 0.184 


<5x>=0.39, <Sy>= 


U. -L / U 


MCC 


0.389999999988E00 


QR 


0.38988E00 


MCS 


0.333E01 




SVD 


0.389999999988E00 


MCS 


0.39002E00 


MCC 


0.377E01 




GS 


0.389999999997E00 


MCC 


0.39003E00 


GS 


0.392E01 




MCS 


0.389999999998E00 


SVD 


0.39003E00 


QR 


0.399E01 




QR 


0.389999999999E00 


GS 


0.39004E00 


SVD 


INF= -1 




TRM 


0.390000429338E00 


TRM 


0.31130E06 


TRM 


0.137E12 




<6x>=0m, <<5y>= 0.320 


<6x>=0.Q0, <(5y>= 0.282 


<,5^>=0.60, <6y>= 


0.269 


MCC 


0.597599999976E00 


MCC 


0.59754E00 


MCS 


0.197E01 




SVD 


0.597599999976E00 


SVD 


0.59754E00 


MCC 


0.334E01 




MCS 


0.597599999983E00 


QR 


0.59761E00 


QR 


0.430E01 




QR 


0.597599999985E00 


GS 


0.59762E00 


GS 


0.471E01 




GS 


0.597599999990E00 


MCS 


0.59763E00 


SVD 


INF= -1 




TRM 


0.597600286570E00 


TRM 


0.19820E06 


TRM 


0.377E11 
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The notation used in Figures 1-9: 



□SM| - 
TRM - 
MCS| - 
QR - 
GS - 
SVD - 
MCC - 




t = -13 -12 -11 -10 -9 -8 -5 -2 -1 3 

Fig.l. At /i(PF)=0.533E07 — being well-posed, Ni = 117. 




t= -7-6 -5 -4 -2 -1 10 12 13 

Fig. 2. At j2{W)=0M5E15 — being ill-posed, N2 = 43. 



TRM - 
QR - 
GS - 
0SM| - 
MCS| - 
MCC - 




t= -2 -1 1 2 3 10 17 20 23 54 73 

Fig. 3. At jl{W) >0.450E16 — beind pathologically ill-posed, N3 = 46. 



□SM| - 
QR - 
GS - 
MCS| - 
SVD - 
MCC - 




t = -14 -13 -12 -10 -8 -7 -1 

Fig. 4. At fi{W)=0.857m6 — being well-posed, N5 = 16. 
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"T 1 1 T" 



t= -7 -6 -5 -3 -2 1 11 

Fig. 5. At /x(VF)=0.308E15 — being ill-posed, Nq = 3. 



10* 




t= -2 -1 1 3 4 17 t = -14 -13-12 -9 -7 -6 1 74 

Fig. 6. At fl{W) >0.450E16 — being patho- Fig. 7. At fi{W) =0.675E07 — being well- 
logically ill-posed, TVy = 5. posed, Ng = 39. 




t= -7-6 -5 -3 -2 74 

Fig. 8._ At fl{W) =0.589E15 — being ill- 
posed, iVio = 2. 




t = -3 -2 -1 1 2 4 74 

Fig. 9. At fl{W) >0.450E16 — being patho- 
logically ill-posed. Nil = 7. 



Fig. 10. 
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0.60- 



0.50- 



0.40- 



0.30- 



0.20- 



0.10- 



0.00- 



: 



0.60- 



0.50- 



0.40- 



0.30- 



0.20- 



0.10- 



<6y> 



0.00 



f I I 



1 



() 



I 



<(5y> 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.00 0.05 0.10 0.15 0.20 0.25 0.30 

Fig.ll (3<m<6, 0.524E03</x(A)<0.150E08) Fig.l2 (7<m<ll, 0.475E09</x(^) <0.518E15) 



5.00- 
4.00^ 
3.00 
2.00 
1.00 



0.00- 



I 
I 

A 

I 
I 
I 

® 

I 
I 
I 

i 

I 
I 

_L_ 



A 
I 

I 

I 

I 

I 

t 

I 
I 

¥ 
I 
I 
I 
I 

_L_ 



I 
I 

I 

? 
I 

I 

I 
I 

# 

I 
I 
I 

_L, 



A 
I 
I 
I 
I 

4 
f 



I 
I 

I 

I 
I 

4 

I 
I 
I 

I <dY> 

_U 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 



The notation: 

X - MCC 
• -MCS 
■ - GS 

o - QR 
□ - SVD 
A - TRM 
la -( □ = x) 



«5^>=< 



l|AX|| 



<Sy>=< 



W 

l|Ay|| 



> 



> 



Fig.l3 (12< m <13, n{A) >0.450E16) 
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The analysis of numerical results reported in Tables 1 -j- 13 and their graphical inter- 
pretation with the use of Figs. 1 13 show that our programs MCC and MCS provide, on 
the average, better accuracy characteristics as compared to the most known analogous 
programs. 

The program MCC has also better time characteristics in the case W — €2, no matter, 
whether a system of equations is well- or ill- or pathologocally ill-posed, but MCC and MCS 

are about twice as worse in time as the program GS (DBEQN) in the case W — C3. This 
is owing to the time consumption on the analysis of zeros in computing 5- ■ — elements 
of matrices B = C3 and on testing various inequalities in accordance with the algorithm 
(2.1) (2.5). The programs MCC and MCS work about 10 times as slow as the program 
GS (DEQN) in the general case W : A — , A . This is due to considerable time 

consumption on reduction of the system WX = y of the general form to systems of the 
type (1.2) and (1.3). 

^From the analysis presented it follows that the critical-component method in its 
qualitative characteristics is the best one of the methods of solution of degenerate and 
ill-posed systems of linear algebraic equations. 

4. Conclusion 

In this paper, we have demonstrated the efficiency of the critical- component method 
for numerical solution of degenerate and ill-posed systems of linear algebraic equations. 

We have proved the theorem according to which the only stable normal solution can 
surely be obtained for degenerate and ill-posed systems of linear algebraic equations by 
the critical-component method. 

Results of numerical experiments (278 examples were computed) on the calculation of 
basic characteristics of solution of the system WX = Y are presented, and a comparative 
analysis has been performed, which shows that the programs MCC and MCS have, on the 
average, better characteristics. 
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5. Appendix 

/. Test examples of systems of equations C2X — Y with two-diagonal matrices of the 
general form: 
System 1 



Co 



1 2 



1 2 
1 



i = l,2,...,M, 



_ 3i+l 

1,2,...,M-1; 



, yj^ = 1/M, 



System 2 



£n r 



£n r 



x, = l/(2f + £*). 2 = 1,2,...,M, 

yi = (2i+£*)(2i+£;+2)' = ^0/ (2Ai' + £o)> 

^ = 1,2,...,M-1, 
where r = 1 - Sq, Sq = Q, 01; 



System 3 



Co 



I il 

5 3 



I ii 
5 3 
7 
5 



x, = l/(2i + l), 
i = l,2,...,M, 

i = l,2,...,M-l; 



System 4 



Co 



5^2 
-1 2 



-1 2 
£t 2 
-1 2 



-1 2 

^1 



a;. = (-l)^+ia, 
2 = 1,2,...,M, 

2/1 = (^0 - 2)a, 

= (-l)^3a, 
i = 2,3,...,A;-l,A; + l,...,M-l, 

2/, = (-l)^(2-£*)a, 

where a^l+e*^, £^ = 0, 0000001, 
el = 0,0001; 



System 5 



Co = 



3 7 



3 7 
3 



a^i = 

z = l,2,...,M, 

% = Vm = 3, 
i = 1,2,...,M- 1. 
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//. Test examples of systems of linear equations C3X — Y with tridiagonal matrices 
of the general form: 



System 6 



System 7 



System 8 



2 -1 
-1 2-1 

-1 2-1 
-1 2 



-1 1 
1-2 1 



1-2 1 
1 a 



1 -1 
-1 1 -1 



■1 1-1 
-1 1 



System 9 



1 r 
p 1 r 



p 1 r 
p 1 



System 10 



1,2,...,M, y, = 2, 
, i = 2,3,...,M-l; 



M-2 
M(M-l) ■ 



Xi = l + (-l)%*, i = l,2,..,M, 

y, = 2el, = {-ly-Hsl, z = 2, 3, M - 1, 



-0' 



where a = £0 = 0, 0000001; 



1 



1 



~ 2i' * ~ 1' 2' Z/l — 4) Z/m — 2M(1-M) ' 



= 1, i = 1,2,...,M, yi = 2-£*, 
yj;^ = 2 + 4, I/, = 3, i = 2,3,..,M-l, 
where p = 1 + r = 1 - e^, £0 = 0, 0000001; 



6 3 
4 6 3 



4 6 3 
4 6 



= 1, i = 1,2,..., Af, 
7/1 = 9,^/^^ = 10, % = 13, z = 2,3,...,M-l. 
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///. Test examples of systems of equations AX — Y with A ^ A^ - filled matrices of 

the general form: 
System 11 



A 



M M-l M-2 ••• 3 2 333 

M-1 M-l M-2 - 3 2 1 

M-2 M-2 M-2 ■■■3 2 1 



3 
2 



3 2 1 
2 2 1 
111 



Vi 



1/i, i = l,2,...,M, 

^ V-fc+l I 333 „, 

k ~^ M ' fM 

k=l 



M 

k=2^ 



'0; 



M 



(M-i + l)Ei+ E 
k=r k=i+i 

i = 2,3,...,M- 1, 
where e*^ = 0,0000001; 



M-fc+l 



System 12 



i = l,2,...,M-l, j = l,2,...^ 
(J'Mi ~ 333, 

^Ml — M+j-l' ~ ^' 



= l/(2i + l), i= 1,2, 
M ^ 
Z/i = E^ 



I — 
Vm 



k=l 

1,2 



M 

E 

fe=2' 



(2fc+l)(i+fe-l) ' 

M-l, 



M, 



1 



,(2fc+l)(i+fe-l) 



+ 111; 



System 13 



A 
j -- 

O'lM 



v%J> «li 



a 



M-j+l ' 



1,2,..., M-l, 

— 1 + ^Q, ^a^;^ 

'^U ~ S'i ~ M-i+1 ' 

i = 2,3,...,M, J =2,3,...,i. 



1 - £*, i 



= 1,2,...,M, 

M-l 

(1-4)(Em^ + 1 + £S), 



fc=l 



yi = (i-£S)(MZ^ 

i = 2,3, ...,M- 1, 
2/m 

= 0,00001; 



M-fc+l 
M 

+ E 

k=l 
M 



1 



M-fe+1 



where en 



System 14 



M M-l 
M-l M-l M-2 



3 2 
2 2 1 
111 



X, = (-l)Vi, i = 1,2,...,M, 

. (-l)'(i-M) x f (-1)^ 

^ fe=i ^ 

1,2,..., M-l, 

M 

yM = E (-i)V^; 
fe=i 



Vi 
i - 



System 15 



^ (%■)! i-j+M' 

i = l,2,...,M, j = l,2,, 



M, 



1,2, 



: 1/i, i 

M ^ 
^^{i-k+M) ' 



1,2, 



M. 



M, 
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IV. Test examples of systems of equations AX — Y with A — - filled matrices of 
the general form: 



System 16 
A = (a^j), a 



Ij S'l M-j+1 ' 

M- 1, 



1 



J = 1,2,..., 



~ ^ji ~ M-i+l' 

z = 2,3,...,M, j = 2,3,...,t, 
x,= {t + l)/i, i = 1,2,...,M, 



Vi 



E 

k=l 
1 



I (M+£*)(M+1) 
k{M-k+l) M ' 



M 

Vm = E 

k=2 



fc+1 
fc 



fc=j+l 

2(M + e*), 



k{M-k+l) ' 



1, 2, M - 1, where el = 0, 0000001; 



System 17{A - the Hilbert matrix) 



1,2,...,M, j = l,2,...,M, 



Vi 



l/i, t = l,2,...,M, 

M 



k{i+k-i) 



1,2, ...,M; 



System 18 



1,2,...,M-1, 



j 

^Mi = "^1^^ = 333, a,ij = a.ji = j^ijzii, 
i = 2,3,...,M, j = 2,3,...,i, 
x, = {-iy/i, i = l,2,...,M, 



„, _ Ml! + (-1)^333 

K=l 

= Ej^ A:(fc+M-1) + 333; 



System 19 



M-l M-l M-2 •••3 2 1 



2 2 2 
Mil 



2 2 1 
111 



x,^l-e*o, i = l,2,...,M, 

_ (i-£;^)(m2-m+2£*) 
yi — 2 ' 

^ ^ (l-e*)(M-0(i+M-3) 

yj^={2M-\){l-e*), 

i = 2, 3, M - 1, where £^ = 0, 0000001; 



System 20(det(A)=0) 



a 


b 


b ■ 


■ b 


a 


b 


a 


b ■ 


■ b 


b 


b 


b 


b ■ 


■ a 


b 


a 


b 


b ■ 


■ b 


a 



(-1)' 

2i+l ' 



yi = bE 



k=l 



I 

M-l 

k=2 

(-1) 
2A;+1 



1,2,...,M, 



(-1) 



1 ~ "'\2M+1 3j ' 



2i+l ' 



2, 3, M — 1, where a = 1 — £, 



0, 



b^l + s 



0, ^0 



1 • 10-11. 
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